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r.  INTRODUCTION 


This  report  deals  with  the  forces  and  moments  exerted  by  a  liquid  payload 
during  spin-up  on  a  spinning  and  nutating  right  circular  cylinder.  The  cur¬ 
rent  work  steins  from  the  original  paper  by  Stewartson,^  improved  by  the 
viscous  correction  of  Wedemeyer.^  It  represents  an  extension  of  the  work 
presented  in  Reference  3,  which  treats  the  particular  case  of  circular  angular 
motion  wuh  constant  nutational  frequency  and  zero  yaw  growth.  Here  we  deter¬ 
mine  the  viscous  shear  contribution  to  the  moment  in  addition  to  that  of  the 
pressure.  The  ad  hoc  model  of  Murphy,"*  described  in  Reference  3,  is  employed 
to  compute  the  spin-up  moment. 

In  tne  analysis,  two  coordinate  systems  are  considered.  The  first  is  an 
inertial  system,  the  unyawed  reference  frame  shown  in  Figure  1,  in  which 
cylindrical  polar  coordinates*  (r,  0,  x  )  and  Cartesian  coordinates  (y,  z,  x) 
are  used.  The  second  coordinate  system,  used  to  describe  the  projectile 

motion,  is  the  y,  z,  x  non-spinning  system  that  has  the  x-axis  along  the 

projectile  axis  of  symmetry;  they  and  z  axes  are  omitted  from  Figure  1  for 

clarity.  The  x  =  0  and  x  =  0  values  are  located  at  the  midplanes  of  the 

unyawed  and  yawed  cylinders,  respectively.  The  x-axis  is  nutating  about  the 
x-axis  with  the  angle  Ki(t);  the  pivot  point  lies  at  the  midplane.  The 

components  of  the  projection  in  the  y,  z  plane  of  a  unit  vector  lying  on 

the  x-axis  are  denoted  by  ny^  and  n^^,  respectively.  It  is  convenient  to  com¬ 

bine  the  two  components  of  yaw  into  a  single  complex  variable 


i^.^-lny^+in^^).  (1.1) 

The  nomenclature  here  is  that  prescribed  in  Reference  5  and  used  in  References 
3,  6,  and  7. 


*  .  .  .  .  .  . 

Pefvmttons  oj  qiiuntzties  are  given  tn  the  LIST  OF  SYMBOLS  Section. 

I.  K.  StewartsoK,  "On  the  Stability  of  a  Spinyiing  Top  Containing  Liquid," 
■hnirnal  af  Fluid  Mechanics,  Vol.  5,  Part  4,  1959. 

P..  F.  ’I.  Weiemeyer,  "Visaoue  Corrections  to  Stewartson’ s  Stability  Criter¬ 
ion,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland, 
BRL  Report  No.  1287 ,  June  1966.  (AD  489687) 

J,  V.  'levbi.r,  "Coy.tributinn  of  Pressure  to  the  Moment  Muring  Spin-Up  on  a 
Uutitio.g  Ligmd-Filled  Cylinder:  Ad  Hoc  Model,"  Ballistic  Research 
Labor  It  jry ,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02566 ,  June  1984. 
(Al  Al‘jJr76' 

4.  ,-7.  Murphy,  "Moment  Induced  by  Liquid  Payload  During  Spin-Up  Without  a 

Critical  layer,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  ARBRL-TR-CPSHl ,  August  1984.  (AD  A1457  16)  (See  also  AIAA  22nd 
Aerospace  Sciences  Meetini,  Reno,  Nevada,  AIAA  Paper  No.  84-0229,  January 
1984.) 


The  projectile  motion  Is  specified  as 


I  =  e’*^  =  Kq 


where 


K]  =  Kq  =  T$t,  f  =  (1  -  ie)  T. 


(1.2) 

(1.3) 


Here  t  is  time,  ^  (rad/s)  is  the  axial  component  of  the  angular  velocity  of 
the  cylinder;  t  and  ex  are  nutational  frequency  and  yaw  growth  rate, 
respectively,  divided  by  1.  Kg  is  the  magnitude  of  the  yaw  at  time  t  = 

0;  (j)j^  is  the  angular  orientation  of  the  x-axis  in  the  y,  z,  x  system  as  shown 

in  Figure  1.  Lengths  and  distances  are  non-dimensional  1  zed  by  a,  the  cross- 
sectional  radius  of  the  cylinder.  The  yaw  grows  when  et  >  0;  however,  in  this 
study  we  apply  the  restriction  e  =  0. 

Computational  results  will  be  presented  in  terms  of  time  histories  of 
non-dimensional  moment  coefficients.  Each  history  is  specified  by  the  three 
parameters  Re,  A  =  c/a,  and  t,  where 


Re 


(1.4) 


is  the  Reynolds  number  and  v  is  the  kinematic  viscosity  of  the  liquid.  The 
reader  is  referred  to  References  3  and  7  for  details  of  the  analyses;  only  the 
relevant  features  will  be  repeated  here. 


II.  FLOW  PROBLEM 


A.  Flow  Variables. 

Here,  as  in  previous  workf’*’^  we  assume  the  yaw  to  be  sufficiently 
small  so  that  a  linearized  analysis  is  applicable,  i.e.,  the  flow  may  be 


5.  C.  H.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous 
Liquid  Payload,"  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground, 
Maryland,  ARBRL-MR-03194 ,  August  1982,  (AD  A1 18676)  (See  also  Jourmal  of 
Guidance,  Control,  and  Dynamics,  Vol,  6,  July-August  1983,  pp,  280-286.) 

6.  N.  Gerber,  R,  Sedney ,  and  J.  M.  Bart  os ,  "Pressure  Moment  on  a  Liquid- 
Filled  Projectile:  Solid  Body  Rotation,"  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02422 ,  October  1982. 

(AD  A120b67) 


7.  N.  Gerber  and  R.  Sedney,  "Moment  on  a  Liquid-Filled  Spinning  and  Nutating 
Projectile:  Solvi  Body  Rotation,"  Ballistic  Research  Laboratory ,  Aberdeen 
Proving  Ground,  Maryland,  ARBRL-TR-02470 ,  February  1983.  (AD  A12S332) 


considered  ‘he  sum  of  a  known  basic  unperturbed  axisymmetric  flow  and  a  3-D 
perturbation  flow.  The  basic  flow  employed  here  is  the  Wedemeyer  spin-up 
model  ^  (see  Eqs.  (2.2)  and  (2.3)  in  Reference  3)  which  is  obtained  by  finite 


difference  calculations.  Thus, 

u  =  iJ  (r,  Jt)  -  Kg  u  (r,  e,  x.  It)  (2.1a) 

V  =  V  (r,  |t)  -  Kg  V  (r.  9,  x,  ^t)  (2.1b) 

w  =  W  (r,  X,  |t)  -  Kj^  w  (r,  0,  x,  ^t)  (2.1c) 

p  =  P  (r,  |t)  -  Kg  p  (r,  0,  X,  4,t).  (2. Id) 


The  u,  V,  w  are  velocity  components  in  the  radial,  azimuthal,  and  axial 
directions,  respectively;  U,  V,  U  are  the  corresponding  velocity  components  of 

the  basic  flow.  The  velocity  components  of  the  perturbed  flow  are  u,  v,  and 

w.  The  quantity  p  is  pressure,  P  is  the  pressure  of  the  basic  flow,  and  p  is 

the  perturbation  pressure.  Velocity  is  non-dimensional i zed  by  a|  and  pressure 

by  pd~ I  ,  where  p  is  the  density  of  the  liquid.  The  angular  motion  prescribed 
by  Eq.  (1.2)  leads  to  boundary  conditions  on  the  perturbed  velocity  components 

at  the  sidewall  and  endwalls  with  a  t,  0  dependence  of  (Tift-0)^ 

B.  Assumptions  of  Model. 

The  present  model  assumes  the  following  form  for  the  perturbed  flow 
vari abl es : 


u  =  Real  [u  (r,  x;  ^t)  exp  {i  (ilt-h)}]  (2. 2d) 

V  =  Real  [v  (r,  x;  ^t)  exp  {i  (ilt-o)}]  (2.2b) 

w  =  Real  [w  (r,  x;  ^t)  exp  {i  (ij-t-o)}]  (2.2c) 

p  =  Real  [g  (r,  x;  |t)  exp  {i  (Tj.t-9)}].  (2. 2d) 


,  K.  H.  W.' ii’iii.’./rr ,  "The  'h:t'teadn  Flow  Within  a  Spinning;  Cylinder,"  Ballistic 
Research  Lahtrratory ,  Aberdeen  Proainj  Ground,  Maryland,  BRL  Report  No. 
r.’.'-d) ,  October  (AD  A4M846)  (See  also  Journal  of  Fluid  Mechanics, 

Vol.  TO,  Part  ,■> ,  1984,  pp.  88.^-399.) 
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The  u,  V,  w,  and  g  are  complex  functions;  the  time,  t,  enters  them  only 

as  a  parameter  due  to  the  "quasi-steady"  assumption  (discussed  at  the  end  of 
Section  II. A  in  Reference  3).  The  resulting  problem  in  u,  y,  w,  and  g  is 

formally  time-independent,  implying  that  there  are  no  transients  in  the  lag 
between  the  motion  of  the  cylinder  and  the  corresponding  response  of  the 
flow.  The  solution  of  Eq.  (2.2)  would  be  most  valid  late  in  the  spin-up 
history  of  the  fluid,  and  its  applicability  would  be  expected  to  decrease  as  t 
decreases. 

The  differential  equations  for  u,  y,  w,  and  g  are  given  by  Eq.  (2.9)  in 

Reference  3,  the  axis  and  sidewall  boundary  conditions  by  Eqs.  (2.10)  and 
(2.11)  in  the  same  report.  V(r)  and  3V/3r  occur  in  these,  but  U  and  W  are 
absent. 

The  following  ad  hoc  endwall  boundary  condition  is  applied  in  the  absence 
of  a  rationally  derived  boundary  condition: 


L(w)  =  (1  +  6c^  3/3x)  w  =  i  (1  -  t)  r  at  x  =  t  A.  (2.3) 


It  is  discussed  briefly  in  the  last  two  paragraphs  of  Chapter  I  of  Reference 
3.  This  condition  will  not  be  satisfied  exactly  by  the  solution  but  only  to 
within  the  accuracy  of  a  least  squares  fit.  The  complex  "thickness"  6Cr  is 
found  as  follows:  ^ 


,-1/2  ,,^1/2 

(1  -  i 

)  (3  -  t)1/2 

-1/2  r,^1/2 

(1  +  i 

)  (1  +  t)^/2 

(2.-1a) 

(2.^0) 

(2.4c) 


with  the  restriction  (3  -  t)  >  0.  The  "effective"  Reynolds  number  used  here 
is  given  by 

Re^  (t)  =  [2  I  V  (r;  |t)  dr]  Re.  (2.5) 

■'o 


For  solid-body  rotation,  with  V  =  r,  the  bracketed  term  is  equal  to  unity,  and 
Re^  =  Re;  while  the  fluid  is  spinning  up,  V  <  r  and  Rer  <  Re. 


c. 


Flow  Solution. 


The  flow  solution  is  expressed  as 


KF 


NJ 


u  =  0  ^  {  r )  X  + 


sin  + 


k  =  1 


j  =  1 


(2.6a) 


KF 


V  =  V  ( r ;  A  + 

0 


k  =  1 


NJ 

sin  X.  X  +  y  d  .  V  .( r)  sin  p.x 
K  ^  J  j  j 

j  =  1 


(2.6b) 


Kl-  NJ 

w  -  w^(r)  -  W|^(r)  cos  Xj^x  -  X  w^(r)  cos  UjX 

j  =  1 


k  =  1 


(2.6c) 


KF 


NJ 


e  =  pQ(r)  X  +  ^  p^(r)  sin  x^x  +  ^  d^  Pj(r)  sin  u^x, 

k  =  1  j  =  I 


(2.6d) 


where  u^,  U|^,  Uj ,  etc.  satisfy  Eq.  (2.23)  of  Reference  3. 


The  X|^‘s  are  solutions  to  the  functional  equation  (see  Eq.  (31)  of 
Reference  6) 

cos  A  +  X|^  6c^  sin  \  =  0,  (2.7) 

so  that  every  term  satisfies  the  homogeneous  endwall  condition 

L(W|^  cos  X|^x)  =0  at  X  =  ±  A,  where  the  operator  L  is  defined  in  Eq.  (2.3). 

At  r  =  1 ,  u^  =  =  w^  =  f),  so  that  the  uj's  are  eigenvalues  of  Eq.  (2.23)  of 

Reference  3  with  replaced  by  Uj  •  (The  normalization  of  the  corresponding 

eigenfunctions  is  described  in  Chapter  III.C  of  Reference  9.)  Thus,  the 
partial  sol  it  Km  t.iken  trom  Fg.  (2.6), 
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-  II  0 


I  r  )  X  ^  I1|^  sin  \^x,  PTC.  , 


sdtistu's  the  sidewdll  conditions  (Eq.  (2.11)  of  KeTerence  3): 


u  ( r'  ^  1 )  -:  u  ^  { r  1)  =  -  i  ( 1  -  t )  x 


(?.9ci) 


v(r-l)-v  (r-l)=-{V(r=l)-i)x 
-  n  r  ' 


(3.db) 


w  ( r  -  1)  "  w  ( r  -  n  ;  1  ( 1  -  r  1  . 


(3.4c) 


\  i  n,i  1  1  V ,  tJi,‘  pcirt  i  d  I  sol  nt  i  on 


u  11  (  r )  X  1-  7  d  .  u  -  ( r) 
o'  ,1  j'  ' 


Sin  jj .  X 


(3.ina) 


(2.10b) 


IS  li'tt  to  satisfy  the  nnn-honKigonoons  endwall  boundary  condition  of  Eq. 
(2.<1.  Tiio  w,^  of  I'q.  (?.inh)  cannot  satisly  Iq.  (2.3)  i  dont  i  ca  1  ly ;  the 

ct>nii>li>\  d|'s  are  therefore  determined  so  as  to  minimize  in  a  least  squares 

sens.'  the  error  over  the  interval  0  <  r  <  1. 

1).  ilperat  ional  IVo<  edure. 

Ihe  current  operational  procedure  ditters  from  the  one  employed  in 
Keterence  1.  Ihe  new  method  was  surpiested  by  C.  II.  Murphy.  The  first  step  is 
to  laliiildfe  the  C()fii(ilex  (Mqenvaluos  and  the  e  i  ijen  f  unct  i  ons  u^,  v.,  w^,  p. 

as  tii'tore  tor  ,|  1 ....  ,  NJ. 


Itie  next  step  is  to  Solve  the  difterential  ei|uations  for  u^^,  v^^,  w^,  p^^ 
(Iq.  (2,23)  of  Keterencp  3,  with  d  anil  ,  1).  The  or'thonorma  1  i  zat  i  on 

k  |i 

jir'oi  e'.s  of  dffeceiuo  0  is  emfiloyed  to  obtain  thrm*  linearly  independent  cnn- 
f  in  loijs  solutions  (‘'j>  1’^,  VJ,  WJ),  where  i  -  1,  2,  3,  by  applying  the 

tollowiMvi  thret'  independi'nt  std.s  of  boundary  cerutifions  at  r  -  1: 


d, .  n  ! , 


V,(M  J. 


w,  (  1  )  0 


(2.11a) 


U^l) 


(1  - 
TT-r 


W3(l)  =  1. 


(2.11c) 


The  third  set,  Eq.  (2.11c),  is  [i  (1  -  x)]  times  the  boundary  values  of  the 
particular  solution  for  solid  body  rotation.  In  providing  a  known  limiting 
case,  these  latter  numbers  prove  to  be  helpful  in  checking  the  program  and 
assessing  the  accuracy  of  the  procedure.  Finally, 


u  =  k,  U,  +  k,,  U„  +  k,  Uo 
0  11  2  2  J  J 


=  k,  V,  +  k„  V„  +  k^ 

0  1  1  2  2  3  i 


(2.12) 


w  =  k,  W,  +  k«  W„  +  k^  W-,  etc, 
0112233’ 


where  k^,  k2. 


and  k3  are  constants  to  be  determined. 


We  apply 


the  same  condition  on  w^(l)  as  on  w(l);  namely,  Eq.  (2.9c); 


w^(l)  =  i  (1  -  t). 


(2.13) 


Then,  by  application  of  Eq.  (2.11), 


k3  =  i  (1  -  t)  (2.14) 

and 

Wo(r)  =  k^  W^(r)  +  k^  W2(r)  +  i  (1  -  x)  W3(r).  (2.15) 


Upon  substitution  from  Eq.  (2.15)  and  Eq.  (2.10b),  the  endwall  ad  hoc 
boundary  condition  at  x  =  A  becomes 


where 


N]  +  2 

R(r)  ij,  (r)  -  ^  e^  W3(r)  =  0, 

j  =  1 


(2.16) 


e-  =  (cos  Uj  A  +  5c^  Uj  sin  A)  dj 


(j  =  1.2, ...NJ) 


®NJ  >  1  '  '^l’ 


®NJ  +  2  ■  *^2 


(2.17a) 

(2.17b) 
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(2.17c) 


+  2  “  "  ^2^'"^ 

li*  (r)  =  i  (1  -  t)  [W3(r)  -  r] 


(2.17d) 

(2.17e) 


The  error,  R(r),  will  not  be  identically  zero.  We  seek  to  optimize  the 
fit  of  the  solution  to  the  boundary  condition  by  determining  a  set  of  ej's 
which  minimizes  the  error  integral 


9  (e 


-NJ  +  2 


)  H 


(2.18) 


The  procedure  for  calculating  the  Cj's  is  presented  in  Appendix  A.  There  are 
now  NJ  +  2  functions  and  constants  with  which  to  minimize  the  error,  as  com¬ 
pared  with  NJ  functions  and  constants  in  the  method  of  Reference  3.  This 
additional  flexibility  is  expected,  generally,  to  produce  an  improvement  in 
accuracy  in  satisfying  the  endwall  boundary  condition.  The  measure  of  the 
relative  error  adopted  here  is  given  by 


Er  =  [g/ 


^  (r)|^  dr]^/^ 


(2.19) 


The  constants  k^,  k2.  and  k3  are  now  evaluated  from  Eqs.  (2.17b)  and  (2.14). 
Finally,  the  differential  equations  (Eq.  (2.23)  of  Reference  3  with 

A  A  A  A 

Cp  =  0)  are  solved  for  U|^,  V|^,  W|^,  p|^  with  the  boundary  conditions  at  r  =  1 


obtained  by 

subtracting  Uq(1)x,  Vg(l)x,  ^1^(1)  from  the  conditions 

of  Eq. 

(2.9).  Thus 

for  k  *  0, 

U|^(l)  =  -  [u^(l)  +  i  (1  -  t)]  b|^ 

(2.20a) 

;^(1)  =  -  Cv'gd)  +  V^(l)  -  x]  b^ 

(2.20b) 

w^(1)  =  0, 

(2.20c) 

where 

Cl  +  i\  «c^)^]  sin  A 

(2.21) 

A  [1  +  (Xj^  6c ^)  3  -  5c £ 
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yiven  by  Eq.  (2.29)  of  Reference  3.  The  complete  solution  is  now  available  to 
compute  forces  and  moments. 


III.  LIQUID  MOMENTS 


A.  Sidewall  Moment. 

Definitions  of  sidewall  pressure  and  shear  moment  coefficients  are  stated 
in  Eq.  (b.6)  of  Reference  7.  Mj^^l  moment  about  the  (transverse)  z  - 
axis  exerted  by  the  liquid  on  the  sidewall  of  the  container: 


A  9 

( 2 upa  c  <()  T  K^)  -  ^  ^{LIM)PL  ^ 


(3.1) 


^(LSM)VL  ^  ^(LIM)VL 


The  LSM  and  LIM  designate  side  and  in-plane  moments,  respectively;  i.e.,  the 
moments  tending  to  change  the  yaw  angle  and  nutational  frequency,  respec¬ 
tively.  P  indicates  pressure,  V  indicates  viscous  wall  shear,  and  the  final  L 
(lateral)  designates  sidewall. 

The  boundary  layer  assumptions  (p.  20  of  Reference  7)  are  retained  in  the 

★  ★  ★ 

derivation  of  the  shear  moment;  thus,  3u/3r,  8u/80,  and  3u/3x  at  r  =  1  are 
dropped.  For  spin-up,  additional  order  of  magnitude  assumptions  are  made. 

3/2 

While  terms  of  0(KQ/Re)  are  kept,  terms  of  0(KQ/Re  )  are  omitted;  thus, 

IJ,  3U/3r,  3W/3r,  and  3W/ 32  are  dropped.  The  term  [Vp(l)  -  1]  appears  in  the 

moment  coefficient  formula,  and  it  is  also  neglected.  Because  of  the  boundary 

layer  in  the  perturbed  flow,  it  is  expected  that  j3v/3rj  and  |3w/3r|  >> 

[3V/3r  -  V/r]  at  r  =  1  except  at  very  early  time  when  the  theory  is  not 
appl ied. 


We  replace  Eq.  (5.7)  of  Reference  7  with  the  following  set  of  formulas 
applicable  to  spin-up  for  the  present  circumstances  (e  =  t  =  0),  where 
integrands  are  evaluated  at  r  =  1: 


^(LSM)PL 


(2  I  A)-^ 


Real 


(g+x) 


dx] 


(3.2a) 


^(LIM)PL  ■ 


A 

Imag  [-i  f  x  (g+x)  dx] 

'-A 


(3.2b) 
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,  A  3W  3v  2 

LSM)VL  "  2  T  A'"Re  ^  IF  "  (-gp  -  t)  }  dx] 

J  -A 


(3.2c) 


.  A  3^  3v  2 

‘  -rmr?  ‘"’^sWf  (■5F-T))dx].  (3.2d) 

J-f 


The  3\//3r  term  replaces  unity  in  Eq.  (5.7c)  of  Reference  7  as  a  consequence  of 

★ 

the  sidewall  condition  on  v  (Eq.  (2.6b)  in  Reference  3).  (The  third  term  of 
the  integral  in  Eqs.  (3.2c)  and  (3. 2d)  has  been  retained  here  and  in  Reference 
7  even  through  other  terms  of  the  same  order  in  Re  have  been  dropped  in  the 
boundary  layer  approximation.)  The  variables  g,  v,  and  w  are  evaluated  by  the 

formulas  in  Eq.  (2.6),  Detailed  expressions  for  the  integrals  in  Eqs.  (3.2a) 
and  (3.2b)  are  given  in  Eqs.  (3.8),  (3,9),  and  (3.10)  of  Reference  3. 

B.  Endwall  Moment. 

Ml2e  is  moment  about  the  (transverse)  z  -axis  exerted  by  the  liquid 
on  the  endwalls  of  the  container.  We  restate  Eq.  (5.12)  of  Reference  7: 


A  0 

4*  "  ^(LSM)PE  *  ^(LIM)PE  ^ 

^(LSM)VE  *  ^(LIM)VE 

The  endwall  pressure  moment  coefficients  are  given  by 

C(lsm)pe  =  ^eal  [i  f  r  {g(r,A)  +  A  V  /r)  dr] 

■^0 

^(LIM)PE  "  +  A  dr], 

'o 

which  are  the  same  as  those  obtained  from  Eqs.  (3.12)  -  (3.18) 


(3.3) 


(3.4a) 


(3.4b) 


in  Reference  3. 


For  the  endwal 1  shear 


^(LSM)VL 


''(LIM)VE 


the  same  as  Eqs.  (5.13c)  and  (5.13d)  of  Reference  7.  The  treatment  of  the 
term  3  (v  -  iu)/3x  parallels  that  in  Section  V  of  Reference  7,  where  it  is 

shown  that  the  second  term  dominates  the  right-hand  side  of  Eq.  (5.19). 
Analogously,  here 


[l/(  T  Re)]  Real 


[l/( T  Re)]  Imag  [ 


Jo 


(3  (v  -  i  u)/3x}  dr] 
X  =  A 


(3.5a) 


{3  (v  -  iu)/3x}  dr],  (3.5b) 
X  =  A 


r  r  (3  (v  -  iu)/3x}  dr  "  Br  (  r  [-A{]-?t  +  V  (r)}  +  i  r  (r)]  dr,  (3.6) 
Jo  ■  ‘  ^  =  A  J  r  V 


where 


KF 


NJ 


k  -  1 


j  =  1 


C.  Total  Moment  Coefficient. 


The  total  side  moment  coefficients  due  to  pressure  and  shear  stress, 
respectively,  are 


^'(LSM)P  ^(LSM)PL  ^(LSM)PE 

(3.8a) 

‘(LSM)V  "  ''•(LSM)VL  ^  C(LSM)VE- 

(3.8b) 

The  total  side  Mioment  coefficient  is 

^'(I  SM)  "  ‘'L5M)P  ^'1LSM)V 

(3.9) 

We  now  lief  me  the  following  side  moment  coefficient: 
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r  =  c 

'  LSM 


T. 


(3.10) 


For  given  cylinder,  liquid,  and  spin,  C|_  is  proportional  to  side  moment, 
regardless  of  the  perturbing  frequency. 

IV.  COMPUTATIONS 


A.  Results. 

We  present  side  moment  coefficient  histories  for  the  same  four  combina¬ 
tions  of  Re  and  A  that  were  considered  in  Reference  3.  These  cover  a  Reynolds 

number  range  of  approximately  5  x  10'^  to  2  x  10^.  Here  we  plot  instead  of 

Clsm"  There  is  an  advantage  to  this  representation  in  that  relative  ampli¬ 
tudes  of  moment  are  immediately  apparent  for  all  perturbing  frequencies. 
Again,  the  results  are  compared  with  inviscid  perturbation  calculations 
produced  by  the  method  of  Murphy.** 

In  Figure  2,  histories  are  shown  for  five  nutational  frequencies: 

r  =  n.04,  0.05,  0.09,  0.12,  and  0.14.  These  were  computed  with  the  Wedemeyer 
laminar  Ekman  layer  spin-up  flow  model.®  The  results  are  qualitatively  simi¬ 
lar  to  those  of  Figure  2  in  Reference  3.  The  x  =  0.04  curve  asymptotically 
approaches  its  maximum  with  increasing  Jt,  indicating  a  resonance  for  spun-up 
liquid.  All  the  other  curves  have  peaks  demonstrating  resonance  with  the  (k, 
n)  =  (3,1)  natural  oscillation  mode  of  the  liquid.*  The  x  =  0.12  and  0.14 
curves  have  additional  peaks  (right-hand  peaks)  produced  by  resonance  with  the 
(k,  n)  =  (6,2)  mode.  There  appears  to  be  a  maximum  overturning  moment  that 
the  casing  can  experience,  produced  by  a  particular  nutational  frequency  at  a 
particular  time.  In  this  instance  the  maxi^mizing  frequency  is  very  likely 
close  to  X  =  0.09  and  the  time  close  to  It  =  1200.  For  x  =  0.12,  the 
amplitude  of  the  (k,  n)  =  (5,2)  mode  is  small  compared  to  that  of  the  (k,  n)  = 
(3,1)  mode;  however,  for  t  =  0.14  the  two  amplitudes  are  comparable. 

The  inviscid  perturbation  results  are  also  shown  for  x  =  0.04,  0.05,  and 
0.09;  for  r  =  0.12  and  0.14  inviscid  results  are  unobtainable  because  of  the 
presence  of  the  "critical  layer"  (discussed  in  References  10  and  11),  occur¬ 
ring  when  V/r  =  t  at  some  r.  The  relative  positions  of  the  curves  are 
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similar  to  .iiose  in  Figure  2  of  Reference  3,  but  the  relative  discrepancies 
are  now  larger,  for  a  discrepancy  in  viscous  shear  moment  is  now  added  to  that 
in  pressure  moment.  However,  the  predictions  of  times  of  peak  moments  still 
agree  to  within  about  0.07  s. 

Figure  3  deals  with  a  case  of  smaller  spin-up  time  (spin-up  time  defined 
in  Chapter  I  of  Reference  12);  therefore,  the  interval  of  interest  applies  to 
smaller  .|.t.  The  Ekman  layer  is  taken  to  be  laminar  in  the  spin-up  flow  model. 
In  Figure  4  the  eigenfrequency  histories  of  five  oscillation  modes  during 
spin-up  show  the  times  when  resonances  with  the  coning  frequency  t  =  0.15 
occur,  i.e.,  when  C[^  =  0.15.  TTie  time  of  occurrence  of  the  peak  obtained  by 

the  present  method,  |t  =  225,  indicates  a  resonance  with  the  (k,  n)  =  (3,1) 

motle.  The  inviscid  perturbation  method  fails  to  show  this  resonance;  possibly 
the  method  deteriorates  with  decreasing  time  shortly  before  it  breaks  down  due 
to  the  critical  layer.  The  inviscid  computation  shows  a  prominent  peak  at 

It  -  250  in  Figure  3  where  the  frequencies  (computed  with  inviscid  perturba¬ 
tions)  of  the  (k,  n)  =  (5,2)  and  (7,3)  modes  are  approximately  equal  to  0.15. 
The  present  calculation  shows  no  indication  of  resonance  at  this  time,  though 
a  jog  appears  in  in  Figure  3  of  Reference  3.* 

The  present  computation  does  not  extend  below  It  -  200  because  of 

problems  (described  in  Chapter  II. C  of  Reference  3)  which  arise  at  "small" 

times.  These  are  manifested  in  the  difficulty  encountered  in  obtaining 
certain  of  the  Mj's.  While  this  difficulty  is  not  necessarily  insuperable, 

the  required  labor  can  at  times  render  the  effort  infeasible. 

C 

In  Figure  5,  moment  coefficient  histories  are  shown  for  Re  =  1.85  x  10 
for  two  aspect  ratios  differing  by  2.55%.  The  turbulent  Ekman  layer  was 

employed  in  the  spin-up  model  here.  The  sensitivity  of  side  moment  to  c/a  is 
borne  out  by  the  two  widely  differing  responses  to  disturbances  of  the  same 
frei|uency.  Both  cases  can  be  plotted  on  the  same  graph  since  the  times  of 
significant  variation  in  the  two  sets  of  curves  differ  markedly  ( Jt  ^  10,000 
for  A  =  4.9/3  and  ]it  '  40,000  for  A  =  5.100). 

We  first  consider  the  c/a  =  5.100  case.  Figure  4  in  Reference  3  demon¬ 
strates  the  agreement  between  the  present  and  inviscid  perturbation  methods 
for  the  pressure  contribution  alone.  Discrepancies  between  the  two  outputs 
are  noticeably  increased  when  shear  moment  is  included,  particularly  near  the 
times  of  occurrence  of  the  peaks.  In  the  present  method,  the  viscous  shear 
increases  the  overturning  moment  by  as  much  as  25%  at  ^t  =  42,000;  whereas, 
the  inviscid  perturbation  method  generally  yields  a  slight  decrease.  However, 
the  two  curves  approach  each  other  as  Jt  increases,  reaching  a  limiting  value 


:  ii  <  in  Lfu  h'cdi’wcyrr  of  .'put 
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of  =  2.‘^b  ^  in"'*.  The  predictions  of  times  of  peak  moments  differ  between 

the  two  methods  :iy  approx  i  mate  I  y  1.6  s  for  C,  tFiaure  4  of  Reference  3) 

artil  by  approximately  I..-  s  for  C^. 

for  c/ci  '  4. d/3  the  ad(iition  of  shear  monent  does  not  affect  the  inviscid 
output  qualitatively,  as  indicated  by  inspection  of  present  Figure  5  and 
Figure  b  of  Reference  3.  However,  it  drastically  changes  the  picture  for  the 
present  m.  tnod;  the  C|_  curve  has  a  minimum  with  negative  value  as  before,  but 

It  also  has  a  peak  immediately  afterward.  The  peaks  for  both  aspect  ratios 
indicate  resonance  with  the  (k,  n)  ^  (b,l)  mode  of  the  liquid. 

3 .  Aci.iira.-y  of  Results. 


ContidenC''  in  the  numerical  result'  rupn’-cs  some  estimate  of  their 
accuracy .  A  thorough-goi  fig  error  stu;ly  for  widt^  ranges  of  physical  parameters 
is  not  feasible;  we  confine  our  attention  tf)  the  cases  treated  in  Section 
IV. A,  There  are  at  least  twelve  parameters  in  the  spin-up  and  perturbed 
Navier-Stokes  calculations  alone  that  ^lontrol  the  accuracy  of  the  numerical 
procedures.  Accuracy  as  used  here  is  quite  different  from  validation  of  the 
tht'ciry.  Thf're  are,  in  fact,  no  experimental  data  to  compare  with  the  present 
theory . 

[n  all  cases  we  computed  spin-up  profiles  with  r-ir.terval  A  r  -  1/2H0  and 
time  interval  A  (it)  small  enough  so  that  five  or  fewer  iterations  were 
requi'“ed  for  convergence.  For  the  perturbation  equations,  we  took  A  r  = 

I  '  lilHi)  and  performed  no  less  than  fifty  orthonormal  i 7at i ons .  We  used  as  many 
a,s  thirty  ti?rins  of  the  power  series  to  evaluate  the  solution  to  the  perturba¬ 
tion  equations  at  r  =  O.UOl  where  the  numerical  integration  began,  and  we  used 
seven  terms  (b|_.  sin  \yx)  in  the  biorthogonal  expansion  of  the  function  x.  In 

all  iterative  proct'sses,  we  required  convergence  at  least  to  five  decimal 

pi  (h  f'S  . 

In  this  proiileri  an  adiiitional  factor,  apart  from  those  introiluced  by  the 
niiiierical  opiirations,  must  be  considered;  namely,  the  error  Er  of  Eq.  (2.19) 
'■■'s.i  1 1 1 'ig  from  approximat  i  ng  ;  (r)  by  a  finite  series,  T.  e.  w^  (r),  in 

,1  =  1  ■’ 

iqs.  (  .’.Ibl  and  This  approximation  is  aiaie  in  arder  to  minimise  the 

V’evill  enr'or  resultincj  from  failure  of  the  Oilurion  to  satisfy  the  endwal! 

'  ounC.iry  conilition  identically.  In  the  pror.edu''c  •i*'  Reference  3,  designated 

“V-m  >',)  I  1"  here,  we  used  Nd  w.'s;  whereas,  in  in.-  |irt-ae''t  procedure,  desig- 
no'ei!  "Metnnd  II,"  we  emj'l  ly  ( \J  ♦  21  ,  :  he  fi-ir.  N'  nt  which  are 

iilentical  to  those  ot  Method  I.  Whenever  possibl.',  we  s^t  N.i  =  6,  but  it  is 
not  always  possible  to  do  so. 

figures  6,  7,  8  and  9  present  compari  sin--,  ot  rorc  for  Methods  !  and  II. 
In  IrOh  methods,  ;  (rl  0  and  all  e ,  '  i)  i*  s-.d:l  bii  :v  rotation,  so  that  in 

jeioi  ipli'  Fr  f  ‘t  ^  ')  IS  i  ndi’t  eraii  nat  Res,r''s  -n-mat.'  a  -jenera’  treno  in 

w"i  ‘”  fdie  error  dei  r'eases  as  tins'  decr.Mses  ‘■’"'r  ‘  r.vichns  a  :’dnim;m',  and 


then  increds  s,  sometimes  at  a  rapid  rate,  as  t  0.  In  Figure  6,  Method  II 
yields  a  considerably  smaller  error  than  Method  I,  being  always  less  than  1/2 
percent;  the  same  result  holds  for  i  =  0.04  and  0.09.  One  would  expect  a  more 
accurate  approximation  from  Method  II  because  there  are  NJ  +  2,  as  compared  to 
NJ,  functions  with  which  to  approximate  !(,■  (r).  This  is  true  in  most  of  our 
cases.  However,  the  error  of  Method  II  can  also  become  very  large,  as  shown 
in  Figure  /.  We  note  that  (r)  and  the  functions  with  which  it  is  fitted 
generally  have  sharp  dips  or  spikes  near  r  =  1  because  of  the  sidewall  bound¬ 
ary  layer.  Furthermore,  at  early  time  when  a  critical  layer  is  present,  these 
functions  often  have  sharp  oscillations  in  the  critical  layer.  These  two 
factors  should  account  in  part  for  the  large  error  of  Method  II  at  early  time. 


Method  II  produces  a  considerably  smaller  error  than  Method  I  for  the 
case  in  Figure  8  and  for  the  c/a  =  b.lOO  case  in  Figure  9.  For  the  c/a  = 

4.9/3  case  in  Figure  9  the  region  of  resonance  is  at  early  time,  and  the  two 
methods  have  comparable  errors.  In  all  the  error  plots  except  those  in  Figure 
1,  six  eigenvalues  and  corresponding  eigenfunctions  were  used.  In  the  plot 

for  T  -  0.12,  due  to  the  difficulty  of  computing  other  eigenvalues,  only  the 
n  =  1,2, 3, 4  modes  were  used  for  {it>1000  and  the  n  =  1,2,4  modes  for  Jt<1000. 

The  eigenvalues  ,i,  together  with  their  associated  eigenfunctions,  w(r), 
are  important  parameters  of  the  solution,  not  only  in  the  number  used  but  also 
in  the  choice  of  eigenvalues  and  functions  when  a  choice  is  available.  At 

solid  body  rotation  the  radial  oscillation  modes  are  identified  according  to 
the  increasing  value  of  the  real  part  of  the  eigenvalue  and  designated  as  n  = 
1,2,  etc.  The  oscillation  modes  retain  their  identity  with  decreasing  time, 

but  eventually  their  history  has  to  be  traced  in  small  time  decrements  since 

some  of  the  Ileal  (  vs  t  curves  cross. 

Lxperience  indicates  that  the  most  accurate  calculations  are  obtained 
when  the  lowest  possible  modes  are  employed  in  the  approximation  to  (r); 
e.g.,  if  Nd  =  4,  th..>  j  =  1,  2,  3,  4  in  Eq.  (2.6)  should  correspond  to  n  =  1, 
3,  4.  Table  1  provides  an  illustration.  To  three  figures,  the  (n  =  1,  2, 
3,  4)  combination  gives  the  same  C|^  values  as  the  six-mode  calculation,  which 

presumably  is  the  most  ac(  urate  choice;  and  the  C|  from  the  (n  =  1,  2,  3) 

combination  differ^  tri'in  tnese  by  less  than  3'.  It  should  be  noted,  however, 
that  agreeiiient  lietwren  tiie  four-  and  si*-modt‘  approximations  is  not  always  as 
good  as  thaf  "Ut.  lined  in  this  particular  case.  The  other  combinations  give 

values  ut  i',!  rn.if  oiffer  by  at  least  Id  and  181  from  0.0204.  In  this  instance 

the  sens  1 ;  1  ^  1 !  y  '  ,  (tioi.e  nt  modes  is  inore  n  itable  than  sensitivity  to  number 
of  mode-. .  m,.  Mill  1  arid  fourth  comi  Illations  in  Table  1  can  arise  in  the 
actual  iimifpr  c  1 1  'i  1  at  i  on.  dnder  certain  circumstances,  the  iteration  process 
may  converg.'  to  an  eigtuivalue  otner  than  tiie  one  In'ing  sought  if  the  initial 
giu'ss  is  n  i!.  .stfuiently'  close  to  the  final  answer.  This  phenomenon  is 
discussed  on  (lage  dl  of  Keterence  3.  In  t  ne  case  shown  in  Table  2,  the  two 
dt'termi  na*  I  or’s  of  j  differ  by  al'ouf  ’•  . 

li'jures  1  :,  il,  and  12  compare  the  ' s  obtained  by  Methods  I  and  II. 

For  the  c’e  and  ke  -  1.8b  ■  Id*’  ;  ases  there  is  general  agreement 

betwc'en  the  outputs  of  the  two  methods  ev(‘n  though  tlie  !  r- ' s  disagree  (Figures 
6  anrl  9).  Ihus,  ir  is  no'  a  sensitive  indicator  ot  the  effect  on  Ci  of  the 


method  employed.  For  the  Re  =  4974  case,  the  discrepancies  are  appreciable; 
e.g.,  at  iit  =  400  the  C[_  ot  Method  I  is  less  than  half  of  the  Cl  of  Method 

II.  In  light  of  the  errors  plotted  in  Figure  8,  we  conclude  that  Method  II 
yields  the  more  accurate  value  of  Cl  in  this  particular  case. 

TABLE  1.  EFFECT  OF  RAOIAL  MODES  ON  ERROR  AND  MOMENT  COEFFICIENT: 

Re  =  4974,  A  =  3.30,  t  =  0.15,  |t  =  220. 


Mode  (n) 
Combi  nations 

Er 

Cl 

,  2,  3,  4,  b. 

6 

0.14 

0.02044 

,  2.  3,  4 

0.58 

0.02043 

,  3,  4,  5 

3.65 

0.02268 

,  3,  4,  6 

10.89 

0.01674 

,  2,  3 

1.38 

0.02051 

TABLE  2.  Er 
A  = 

AND  Cl 
=  3.12, 

FOR  TWO  SETS  OF  MODES: 

T  =  0.14,  $t  =  1000. 

Re  =  39,772, 

Mode  (n) 
Combinations 

Er 

Cl 

.  2,  3,  4,  5, 

6 

0.97 

0.000637 

.  2,  3,  4 

4.14 

0.000585 

V.  SUMMARY 

The  primary  mechanism  for  producing  unstable  motion  of  projectiles  with 
fully  spun-up  liquid  payload  is  resonance  between  the  coning  motion  of  the 
projectile  and  inertial  oscillations  of  the  contained  liquid.  Resonance  also 
occur'',  early  in  the  flight  during  spin-up,  but  the  transient  natures  of  the 
fluid  flow  and  projectile  angular  motion  make  prediction  of  instability  more 
difficult  than  for  solid  body  rotation.  This  report  describes  a  procedure, 
based  on  linearized  flow,  for  calculating  the  side  moment  exerted  by  the 
liquid  on  the  casing  during  spin-up  for  the  special  case  of  constant  coning 
frequency  and  zero  yaw  growth.  Flere  the  moment  due  to  shear  forces  is 
included  in  addition  to  that  due  to  pressure  forces.  The  heuristic  approach 
of  Murphy'*  is  used  which  makes  crucial  assumptions  regarding  time-dependence 


of  the  flow  and  the  endwall  boundary  condit,.:i.  Because  of  the  specialized 
angular  motion  the  present  treatment  cannot  simulate  actual  flight,  although 
it  can  simulate  perforniable  gyroscope  experiments.  The  validity  of  this 
method  is  expected  to  be  strongest  at  late  time  and  diminish  with  decreasing 
time. 

Calculations  were  performed  for  four  combinations  of  Re  and  A  and  a 
variety  of  nutational  frequencies.  The  peaks  of  the  vs  4it  curves  indicate 

transient  coincidence  of  certain  inertial  frequencies  of  the  liquid  with  the 
frequency  of  coning.  For  a  given  container,  payload,  and  spin,  there  is  a 
finite  frequency  which  maximizes  the  side  moment  during  the  course  of  spin-up. 
The  duration  of  a  side  moment  surge,  as  well  as  its  amplitude,  will  undoubted¬ 
ly  be  a  factor  in  determining  onset  of  instability  in  a  free-flight  situation. 
Comparison  of  output  with  results  of  the  inviscid  perturbation  method  shows 
qualitative  agreement  at  late  time  but  differing  indications  of  resonance  at 
early  time,  as,  for  example,  in  Figure  3. 

In  the  current  analysis  an  ad  hoc  correct^  m  was  made  to  the  endwall 
boundary  condition  which  is  an  extension  of  the  correction  for  the  solid  body 
rotation  problem.  A  subsequent  refinement  should  be  the  inclusion  of  the 
Ekman  layers  in  the  basic  flow.  This  step  is  necessary  in  order  to  have  a 
rational  approximation  to  the  solution.  It  is  not  expected  to  be  easy  to 
achieve. 

ACKNOWLEDGEMENTS 

The  author  is  pleased  to  acknowledge  the  help  of  Ms.  Joan  M.  Bartos,  who 
programmed  and  performed  the  moment  calculations.  Appreciation  is  extended  to 
Ur.  Raymond  Sedney  for  his  overall  guidance.  The  moment  data  for  inviscid 
perturbation  were  obtained  from  calculations  performed  with  a  program  prepared 
by  Mr.  James  W.  Bradley. 


Figure  1.  Diagrams  of  Coordinates  and  Cylinder 


ive  Njtational 
734  rad/s. 


Present  Calc. 


Hornent  Coefficient  Mistori 
1.85  X  10^,  ^=541  rad/s,  t 


rison  Between  Percentage  Errors  of  ilethods 
39772,  A  =  3.12,  t  =  0.05. 


gure  7.  Comparison  Between  Percentage  Errors  of  f.ethods  I  and  II 
Re  =  39772,  A  =  3.12,  x  =  0.12. 


O  Method  I 


Between  Percentage  Errors  of  Methods 


Comparison  Bet 
Re  =  4974,  A  - 


son  Between  Side  Moment  Coefficients  of  Methods 


REFEKENCES 


1.  K.  StfWdrtsoii ,  "On  the  Stability  of  a  Spinning  Top  Containing  Liquid," 
Jnurnal  of  Fluid  Mechanics,  Vol.  6,  Part  4,  iqb'J. 

Z,  £.  H.  Wedeineyer,  "Viscous  Corrections  to  Stewartson's  Stability  Criter¬ 
ion,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland, 
BRL  Report  No.  1287,  June  1J66.  (AD  489637) 

3.  N.  Gerber,  "Contribution  of  Pressure  to  the  Moment  During  Spin-Up  on  a 
Nutating  L i qui d-Fi 1 1 ed  Cylinder;  Ad  Hoc  Model,"  Ballistic  Research 
L^a^oi^^^c^i^^^^Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02563,  June  1984. 

4.  C.  H.  Murphy,  "Moment  Induced  by  Liquid  Payload  During  Spin-Up  Without  a 
Critical  Layer,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  ARBRL- TR-025S1 ,  August  1984.  (AD  A145716).  (See  also  AIAA 
?2nd  Aerospace  Sciences  Meeting,  Reno,  Nevada,  AIAA  Paper  No.  84-0229, 
January  1984.) 

b.  C.  H.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous 
Liquid  Payload,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  ARBRL-MR-03194 ,  August  1982.  (AD  A118676)  (See  also  Journal  of 
Guiilance,  Control,  and  Dynamics,  Vol.  6,  July-August  1983,  pp.  280-286. ) 

h.  N.  Gt'rber,  R.  Sedney,  and  J.  M.  Bartos,  "Pressure  Moment  on  a  Liquid- 
Mi  led  Projectile:  Solid  Body  Rotation,"  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  Maryland,  ARBRL- TR-02422 ,  October  1982. 

(AD  A 12036 7) 

N.  (lerber  and  R.  Sedney,  "Moment  on  a  Liquid-Filled  Spinning  and  Nutating 
Projectile:  Solid  Body  Rotation,"  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  ARBRL- TR-02470,  February  1983.  (AD  A125332) 

8.  E.  H.  Wedemeyer,  "The  Unsteady  Flow  Within  a  Spinning  Cylinder," 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  BRL 
Report  No.  1226,  October  1963.  (AD  A431846)  (See  also  Journal  of  Fluid 
Mechanics ,  Vol.  20,  Part  3,  1964,  pp.  383-399.) 

1.  C.  W.  Kitchens,  Jr.,  N.  Gerber,  and  R.  Sedney,  "Oscillations  of  a  Liquid 

in  a  Rotating  Cylinder:  Part  I.  Solid-Body  Rotation,"  Ballistic  Research 

lahoratO''y,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02CB1 ,  June  1978. 

(All  AOS  7/69) 

!.  R.  Sedney  and  N.  Gerber,  "Oscillations  of  a  Liquid  in  a  Rotating 

Cylinder;  IMrt  11.  Spin-Up,"  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Grounit,  Maryland,  ARBRL- TR-n2489 ,  May  1983,  'AD  A129094) 

1.  k.  .u'itney  and  N.  Gerber,  "A  Study  of  the  Critical  Layer  in  a  Rotating 

Liquid  Payload,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  ARBRL-TR-02682,  August  1984,  (See  also  AIAA  Paper  No.  84-0342, 
January  1984.) 

2.  R.  Sedney  and  N.  iGf^rber,  "Viscous  Effects  in  tbe  Wedemeyer  Model  of  Spin- 
Up  triom  Rest,"  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  ARBR L- TR-02493 ,  June  1983.  (AO  A129SJ6) 


APPENDIX  A 
CALCULATION  OF  e j ' s 


)/ 


APPtNi)IX  A.  CALCULATION  OF  ' s 


Lq.  (^.Ib)  >s  dtp'l  : 

NJ  +■  r' 

'<[r)  =  ^  ;  r )  -  ^  e  .  w  ^  (  r  . . 

J  =  1  ' 

We  now  Oefme  tbe  followinij  elements; 


.  1 


i'  ' 


bunj  [w^  (r)]  Jr 


n  =  1 . NJ  +  2 


ui.j 


J  0 


^  »  conj  [w^]  dr 


(A.l) 


(A. 2) 


n,j=l . NJ+2,  (A, 3) 


where  conj  [w^]  is  the  complex  conjugate  of  w^. 


The  ej's  which  minimize  the  error  integral  g  of  Eq.  (2.18)  are  found  by 
solving  the  following  system  of  linear  complex  equations: 


NJ  +  2 

L 

J  -  1 


■  ^n 


n  =  1 , . . . ,  NJ  +  2 . 


(A. 4) 
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d  cross-sectional  radius  of  cylinder.  Figure  1  [cm] 

A  -  c/a,  aspect  ratio  of  cylinder 

b|^  biorthogonal  coefficients  for  series  expansion  of  x  in  interval 

-A  <  X  <  A,  Eq.  (2.21) 

c  half-height  of  cylinder  [cm] 


LIM)PE 

^(LIM)PL 

'-(LIM)VL 

^(LIM)VL 

SM 

LSM)P 

‘•(LSM)PE 

LSM)PL 

‘'(LSM)V 

LSH)VE 

I.SM)VL 


SM  ^  side  moment  coefficient 

endwall  pressure  in-plane  moment  coefficient,  Eqs.  (3.3) 
and  (3.4h) 

sidewall  pressure  in-plane  moment  coefficient,  Eqs.  (3.1) 
and  (3.2b) 

endwall  viscous  shear  in-plane  moment  coefficient,  Eqs.  (3.3) 
afid  (3.bb) 

sidewall  viscous  shear  in-plane  moment  coefficient,  Eqs.  (3.1) 
and  (3. 2d) 

total  side  moment  coefficient,  Eq.  (3.9) 

pressure  side  moment  coefficient,  Eq.  (3.Ba) 

endwall  pressure  side  moment  coefficient,  Eqs.  (3.3)  and  (3.4a) 

sidewall  pressure  side  moment  coefficient,  Eqs.  (3.1)  and  (3.2a) 

viscous  shear  side  moment  coefficient,  Eq.  (3.8b) 

endwall  viscous  shear  side  moment  coefficient,  Eqs.  (3.3) 
and  (3.5d) 

sidewall  viscous  shear  side  moment  coefficient,  Eqs.  (3.1) 
and  (3.2c) 

natural  oscillation  frequency  of  rotating  liquid/l 
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coefficient  in  sei'ies  given  by  the  3rd  terms  on  right-hand 
sides  of  Lqs.  {3.6a)  -  (2.6d) 

coefficient  defined  in  Eq.  (2.17a) 

relative  error  measure  defined  in  Eq.  (2.19) 

c  (1  -  ie)r,  complex  representation  of  angular  motion,  Eq.  (1.3 

error  integral  defined  in  Eq.  (2.18) 

index  of  axial  eigenfunction  and  eigenvalue 

coefficients  of  combination  of  linearly  independent  solutions, 
(Eq.  (2.12) 

index  of  final  term  of  sin  Xj^x  and  cos  X^^x  in  Eq.  (2.6) 

yaw  amplitude  at  time  t  =  0 

Ky  yaw  amplitude  at  time  t,  Eq.  (1.3) 

nondimensi onal  x  (and  x)  coordinate  of  pivot  point 

sidewall  and  endwal 1  contributions,  respectively,  to  liquid 
moment  about  z-axis  [g  cm^Vs'^J 

index  of  radial  mode  for  ei ()enf requency , 

components  in  the  y,  z  plane  of  a  unit  vector  lying  on 
the  x-axis 


number  of  terms  in  the  series  given  by  the  3rd  terms  of  the 
right-hand  sides  of  E(is.  (2.6) 
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2  *2 

pressure/( pd  I  ) 

p  p 

r,  X  variation  of  perturbation  pressure/ ) 

jth  radial  eigenfunction  in  series  given  by  the  3rd 
term  in  Eq.  (2.6d) 


first  term  in  solution  to  g,  Eq.  (2.6d) 


coefficient  of  sin  X|^x  in  series  given  by  the 
2nd  term  of  the  right-hand  side  of  Eq.  (2.6d) 

p  p 

perturbation  pressure/CK^pa  |  ),  Eq.  (2. Id) 


axi symmetric  unperturbed  pressure/ (pa 
spin-up  model,  Eq.  (2. Id) 


) ,  Wedemeyer 


(1/a)  X  radial  coordinate  in  inertial  coordinate  system 


error  function  defined  in  Eq.  (2.16) 

2  ^ 

Reynolds  number  =  a‘^|/v 


"effective"  Reynolds  number,  Eq.  (2.5) 
time  [s] 

(1/Ca|])  X  radial,  azimuthal,  axial  velocity  components, 
respectively,  in  inertial  frame 

(l/[Kga|])  X  r,  X  variation  of  perturbation  radial, 
azimuthal,  axial  velocity  components,  Eq.  (2.2) 

coefficients  of  sin  p-x  and  cos  p-x  in  series  expansions 

J  J 

for  u,  V,  w,  p  (3rd  terms  of  right-hand  sides  of  Eqs.  (2.6)) 
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u 


k  ’ 


u 


0  ’ 


★  ★  ★ 

U,  V,  w 


U,  V.  W 


U. 

1 


V. 


w. 


X,  y,  z 


X,  y,  z 


coefficients  of  sin  Xj^x  and  cos  X|^x  in  series  expansions  for 
u,  V,  w  (2nd  terms  of  right-hand  sides  of  Eqs.  (2.6)) 

1st  terms  in  right-hand  sides  of  Eqs.  (2.6a),  (2.6b),  and  (2.6c) 

(l/[KQa^])  X  radial,  azimuthal,  axial  perturbation  velocity 
components  in  inertial  system,  Eqs.  (2.1)  and  (2.2) 

(l/[aj])  X  radial,  azimuthal,  axial  velocity  components  of 
axisymmetric  unperturbed  flow,  Eq.  (2.1) 

(i  =  1,  2,  3)  linearly  independent  solutions  which  are 
combined  to  obtain  u^,  v^,  w^,  Eq.  (2.12) 

rectangular  coordinates  in  inertial  system  (x-axis  along 
trajectory)  [length/a],  y  =  r  cos  e,  z  =  r  sin  e 


rectangular  coordinates  in  aerobal 1 istic  system  (x-axis 
along  cylinder  axis)  [length/a] 


Be 


iSC 


E 


'P 

0 


quantity  defined  in  Eq.  (2.4b),  used  in  evaluating 

correction  term  in  ad  hoc  endwall  boundary  condition, 

Eqs.  (2.3)  and  (2.4c) 

=  (1,^t)  X  yaw  growth  per  radian  of  nutation 
=  0  for  X|^  #  0,  =  I  for  \  =  0  (Eqs.  (2.23)  of  Reference  1) 
polar  angle  (azimuthal  coordinate)  in  inertial  system 
eigenvalue  in  axial  problem,  Eqs.  (2.6)  and  (2.7) 
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eigenvalues  occurring  in  series  given  by  the  3rd  terms  of 
the  right-hand  sides  of  Eq.  (2.6);  see  paragraph  preceding 
Eq.  (2.8) 

kinematic  viscosity  of  liquid  [cm^/s] 

vector  describing  angular  motion  of  cylinder,  Eqs.  (1.1) 
and  (1.2) 

density  of  liquid  [g/cm^] 

nutational  frequency  of  cylinder/4. 

spin  rate  of  cylinder  [rad/s],  taken  to  be  positive 

function  defined  in  Eq.  (2.17e),  used  in  Eq.  (2.16) 
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